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Introduction 


This report has two sections that address aspects of the research completed to date. 
The first section presents proof for the existence of a single Poisson's ratio function in 
isotropic linear viscoelastic materials. Besides the proof the first section presents an in 
depth discussion of three dimensional viscoelastic material properties and their relationships 
for linear isotropic and orthotropic viscoelastic materials. The second section presents a 
discussion of the alternate invariant definition as used by Abaqus and as relates to the form 
used by Dr. S. Peng. 

NASA's support and encouragement of this research is greatly appreciated. 
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Poisson's Ratio Function in Linear Viscoelasticity 


2.1. Elements of Three Dimensional Elasticity 

For completeness the discussion must begin at the most fundamental level. To be 
consistent with the development of linear viscoelastic theory it is first necessary to assume 
that displacements are small. Therefore, in the consideration of equilibrium of an 
infinitesimal element it is sufficient to treat the equilibrium of the element in the 
undeformed configuration. From the balance of angular momentum for the infinitesimal 
element we find that the stress tensor must be symmetric, i.e. a- = a }i . Thus, there are 
only six unique components of the stress tensor at any point of a body. 

The next consideration that must be addressed is the measure of the deformation. 
Infinitesimal strain is defined by the expression: 
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dU; 3u^ 
dX; + dx ; 


( 1 ) 


where e tj is the infinitesimal strain tensor and u ; is the displacement in the x ; coordinate 
direction. Following the fundamentals of calculus, = e-; the strain tensor is also 
symmetric with six independent components. 

Following the approach developed in the theory of elasticity the stress is related to 
the strain by the constitutive relationship. In elasticity this relationship is: 

ff ij = Cyu (2) 

where is the stiffness tensor and summation over the repeated indices is implicit as is 
common in tensor notation. The stiffness tensor, C^, is a tensor of rank four and 
therefore represents 81 individual constants. The inverse or compliance form of the 
constitutive equation is: 

£ kl = Sklmn Omn, (3) 

where is the compliance tensor. As with the stiffness tensor, the compliance tensor is 
a tensor of rank four with 81 individual components. The compliance tensor is the inverse 
of the stiffness tensor; therefore: 


Qjld S 


klmn 


where 5 im is the identity tensor. Therefore, there are a maximum of 81 independent elastic 
material properties. Symmetry of the stress tensor and symmetry of the strain tensor leads 
to the requirements that: C ijld = C jikl = C ijIk and = S lkmn = S^, which reduces the 
number of independent components to 36. Writing the stored elastic strain energy, U, as: 
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and making the substitution for the stress tensor from equation (2) yields: 

U = |e ij C ijkl e kl . (4) 

Noting then from equation (4) that there is no loss in generality by assuming the additional 
symmetry of the stiffness tensor such that: C ijkl = C^. Therefore, the stiffness tensor (and 
hence the compliance tensor) has 21 independent elastic material constants. 

The constitutive relations given in equations (2) and (3) with the symmetry 
conditions discussed are general relations with no assumptions on symmetries of material 
behavior. A material in which there are no symmetries of behavior is defined as an 
anisotropic material. Introduction of various levels of material symmetries leads to further 
simplifications of the constitutive relation and reductions in the number of independent 
material constants. Orthotropic behavior or orthorhombic symmetry leads to nine 
independent elastic material constants. Isotropic behavior has two independent constants. 


2.2. Elements of One Dimensional Linear Viscoelasticity 

The foregoing discussion is primarily related to elastic behavior. For the analysis 
of viscoelastic behavior it is necessary to return to equation (2) and recast the relationship 
in that appropriate to viscoelastic behavior. A number of forms for the fundamental 
relation between stress and strain may be defined for viscoelastic behavior. The form 
assumed here is that of the historical convolution integral, i.e.: 

t 

7(0 = J^&J(t-r)dr, (5) 

-QO 


or 


<*t) = 




dr 


G(t-r)dr, 


-00 


( 6 ) 


where y(r) is the strain at the time r, o(j) is the stress at the time r, J(r) is the creep 
compliance function, G(r) is the stress relaxation modulus, t is the current time, and r is 
some time on the interval from -oo to the current time, t. If the stress or strain histories 
are known to be discontinuous at some point in time it is possible to arbitrarily define that 
time to be the reference time, i.e. r = 0. Further, if either the stress or strain history is 



defined to be zero for all times prior to the time r = 0, the historical integrals above 
reduce to the form: 


t 


7 (t) - <7(0)J(t) + ^rlj(t-r)dr, 

(7) 

0 


and 


t 

®<t) - 7(0)G(t) + ^ G(t-r)dr. 

(8) 

0 


At this point it is useful to introduce the notation of * to 
operation. Therefore, equations (7) and (8) may be written as: 

indicate the convolution 

7 = J * dcr 

(9) 

and 


a = G * d7, 

(10) 


respectively. Note that the order is A * dB where the "d" is associated with the time 
derivative of the function B. 

Without proof the following identities are introduced for the convolution operator. 
The proofs are easily obtained. 

Commutativity: Q * dP = P * dQ (I) 

Distributivity: (Qj + QJ * dP = Q, * dP + Q 2 * dP, (II) 


Q * d(P t + ? 2 ) = Q * dPj + Q * dP 2 , (in) 
Associativity: (Q * dP) * dR = Q * d(P * dR), (IV) 

Identity Element: Q * dl = 1 * dQ = Q, (V) 

Zero Element: Q * dP = 0 then Q(t) = 0 or P(t) = 0, V t ,(VI) 

Inverse Element: Q * dQ -1 = 1 (VII) 

Inversion of Products: (Q * dP) _I = P -1 * dQ -1 (VHI) 
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where 1 is the Heaviside step function, i.e. l(t) = {0, -oo <t<0;l,0<t< oo}, 0 is 
the zero function, i.e. 0(t) = 0 V t , and ( -I ) indicates the inverse convolution operation. 

As a consequence of identity (I) the stress strain relations in equations (9) and (10) can be 
rewritten as: 

7 = o * dJ ( 11 ) 


and 


a = 7 * dG. (12) 

These relations with the identities are used in the development of the three dimensional 
constitutive equations for linear viscoelastic materials. 


2.3. Elements of Three Dimensional Linear Viscoelasticity 

The equations presented thus far for the analysis of viscoelastic response of a 
material have been presented in a one-dimensional form. For three dimensional 
viscoelastic analysis it is necessary that these equations be generalized to three dimensional 
behavior. Therefore, generalizing equation (2) for viscoelastic behavior yields: 

a ij = Ciju * d £ u> (13) 

note that in this equation the components are all time dependent functions. Further 
recognize that, as a result of the convolution identities (II) and (III), this equation 
represents the summation of nine individual convolution integrals. The compliance form 
of the viscoelastic constitutive relation is: 


e ld ~ Sklmn * dtf^,. (14) 

Following the discussion presented for the elasticity development, the tensor functions 
C p (t) and Skirrm (0 must maintain the same symmetries relative to the index pairs i,j, k,l, 
and m,n, therefore, C ijkl (t) = C jikl (t) = C ijlk (t) and S ldmn (t) = S lkmn (t) = S klnm (t), V t . Then, 

analogous to elasticity, there are only 36 independent material property functions. For a 
viscoelastic material the internal energy can be written as: 


U = t- e • 
2 1J 




then from (13) 


U = • d 


Su * d£ “ 


The analogy with elasticity is almost complete; by use of convolution identities (I) and (IV) 
the order of the strain tensors may be exchanged in the above expression; hence, there is 
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no loss in generality as a result of the assumption that C ; j U (t) = C^t) and S^^O) = 

Sn^CO, V t . 

To complete the analogy with elasticity it is necessary to determine the relationship 
between C ijk] (t) and S^^t). This relationship is found by substitution of (13) into (14). 
This yields: 


®ij ^-"ijld 


* d 


S u™. * 


and from convolution identities (I) and (TV) this becomes: 


ff y = % * d5 im 5 jn = C. jkl * d 


o_ n * dS 

mn klmn 


^ijkl * 


* dS 


klmn 


= \ a * dC 1 * dS = (j * d J C * dS L . 

1 mn ijkl f klmn ™ i ijkl klmn f 


Therefore: 


ik 


0 = a * dJ 8- 5-„ - C * dS 

mn imjn ,jy klmn 


} 


Finally, from convolution identity (VI) for any general stress history that satisfies the 
conditions discussed, i.e. zero to time t = 0, we find: 


Cjjki * — 5 im 5j n , (15) 

which is similar to the result obtained for elasticity, however, not nearly as simple. 

Symmetries of response that lead to the reduction of the number of independent 
material properties are by definition time independent. Therefore, the results obtained for 
elastic materials are directly applicable to viscoelastic materials. Hence, an orthotropic 
viscoelastic material has nine independent material property functions and an isotropic 
viscoelastic material has two independent material property functions. 

For the investigation of the Poisson effect in viscoelastic materials it will be 
assumed that the material in question is at most orthotropic; later discussions will be 
restricted to isotropic materials. To simplify the discussion the compact notation 
introduced by Tsai for composite materials will be used. Compact notation is defined as 
follows: 
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*1 



r -n 

*11 





*2 



*22 


fi 2 


£22 

*3 

> 

= < 

*33 

► and •* 

«3 

> = « 

£33 

*4 



*23 

£4 


£23 

*5 



*13 


£5 


£,3 

. *6. 



k. ^12 „ 


. fi 6 , 


. C 12. 


( 16 ) 


Then for an orthotropic material in which the axes of symmetry or principal material 
directions are alined with the coordinate directions, i.e. a specially orthotropic material, 
the compact notation for the stiffness and compliance is: 


" C n 

C12 

C,3 

0 

0 

0 " 


Cllll 

Cl 122 

C1133 

0 

0 

0 

C12 

C 2 2 

C-23 

0 

0 

0 


Cl 122 

C 2 222 

C 22 33 

0 

0 

0 

C13 

C23 

C33 

0 

0 

0 


^1 133 

C 2233 

C3333 

0 

0 

0 

0 

0 

0 

C 44 

0 

0 


0 

0 

0 

c 4444 

0 

0 

0 

0 

0 

0 

C55 

0 


0 

0 

0 

0 

C5555 

0 

0 

0 

0 

0 

0 

^66 _ 


0 

0 

0 

0 

0 

^6666 _ 


and 


~ Sn 

6*12 

6*13 

0 

0 

0 “ 


6*1111 

^1 122 

6*1133 

0 

0 

0 

S„ 

6*22 

S23 

0 

0 

0 


^ 1 122 

6*2222 

S 22 33 

0 

0 

0 

$13 

6*23 

6*33 

0 

0 

0 


^ 1 133 

6*2233 

6*3333 

0 

0 

0 

0 

0 

0 

S44 

0 

0 


0 

0 

0 

S4444 

0 

0 

0 

0 

0 

0 

6*55 

0 


0 

0 

0 

0 

6*5555 

0 

0 

0 

0 

0 

0 

6*66 _ 


0 

0 

0 

0 

0 

6*6666 _ 


In contracted notation the interrelation of stiffness and compliance, equation (15), 
becomes: 


C S * dS jt = 5; 


ik’ 


(19) 


where the repeated index indicates a summation from 1 to 6. Expanded in matrix form the 
above is: 
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Cl, 

Cj 2 

C13 

0 

0 

0 

C,2 

C 22 

C23 

0 

0 

0 

C, 3 

C23 

C33 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0 

C55 

0 

0 

0 

0 

0 

0 

^66 


^11 

S 12 

S 13 

0 

0 

0 

■*12 

s 22 

S 23 

0 

0 

0 

5 13 

S 23 

S 33 

0 

0 

0 

0 

0 

0 

S 44 

0 

0 

0 

0 

0 

0 

S 55 

0 

0 

0 

0 

0 

0 

S 6 6 


1 0 0 0 0 0 

0 1 0 0 0 0 

0 0 1 0 0 0 

0 0 0 1 0 0 

0 0 0 0 1 0 

_0 0 0 0 0 1 _ 


( 20 ) 


which represents twelve nontrivial equations interrelating the compliance functions and 
stiffness functions. Recognizing that the matrices may be partitioned into four 3x3 sub- 
matrices, i.e. 

(^Tension q 
0 C Shear 

where C Tens ^ on and are the sub-matrices containing the tensile stiffness functions 

and shear stiffness functions, respectively, S^ ens ^ on and S^hear are the corresponding 
compliance sub-matrices and I is the identity matrix. Written in this form it is clear that 
the off diagonal sub-matrices make no contribution. Therefore, the relations may be 
reduced to two sets of matrix equations, which are: 

^.Tension 
and 

.Shear 


Noting that and S^hear are diagonal matrices, equation (23) reduces to three 

independent equations, which are: 

1 = C* * dS*. (24) 


* d 


..Shear 


= I 


(23) 



* d 


..Tension 


ion _ i 


( 22 ) 


* d 


^Tension q 
0 sShear 


I o' 

.0 ij ’ 


( 21 ) 


1 = C 55 

* dS 55 , 

(25) 

1 = C 66 

* dS 66 . 

(26) 
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Therefore, for a specially orthotropic viscoelastic material the shear stiffness functions and 
shear compliance functions are inverse convolution functions of one another. However, 
the same comment does not hold for the tensile terms. 


2.4. Engineering Properties in Linear Viscoelasticity 

If one defines the shear stress relaxation modulus functions such that: 

= C44 * de 4 = G23 * d7 4 , (27) 

C5 = C55 * de 5 = G 13 * d7 5 , (28) 

a f> = ^66 * de 6 = G 12 * d7 6 , (29) 

where Gy are the shear stress relaxation modulus functions in the i-j plane and 7 k is the 

engineering shear strain in the associated plane, the relationship between the stiffnesses and 
moduli is a one to one relationship. Similarly, by defining the shear creep compliance 
functions such that 


£4 = S44 * d<r 4 = * d tr 4 , (30) 

£5 = S 55 * d(r 5 = J 13 * dcr 5 , (31) 

£ 6 = ^66 * = ^12 * (32) 

an analogous one to one relationship exists between the two sets of compliance functions. 
Then as a consequence of equations (24) - (26) and the above equations (equations (27) - 
(32)), the engineering creep compliance functions and stress relaxation functions for shear 
are also inverse convolution functions, i.e.: 


II 

* dG 12 , 

(33) 

II 

U> 

* dG 13 , 

(34) 

1=123 

* dG23 . 

(35) 


This result is in full agreement with the one dimensional approach to viscoelasticity and the 
results stated by Staverman and Schwarzl 1 for shear behavior. 

Before continuing with the investigation of Poisson's effect in a viscoelastic 
material a discussion of the interrelation of the stiffness and compliance functions is 
appropriate. To begin the discussion it is first important to note that terms Cy in equation 
(22) are not the tensile stress relaxation modulus functions, they are the tensile stiffness 
functions. They are dependent on the tensile (engineering) stress relaxation modulus 
functions as well as other functions which describe the Poisson effect of the material. The 
compliance terms S H are directly the tensile engineering creep compliance functions. 
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Therefore, by defining that F(t) is a tensile creep compliance function, after the notation of 
Staverman and Schwarzl, the creep compliance functions and the tensile components of the 
compliance tensor are: 

S fi (t) = F,(t), (36) 

where F ; (t) is the engineering creep compliance function in the ith principal material 
direction. 

For clarity, consider momentarily a orthotropic elastic material. In an elastic S ;i = 
1/E„ where E, is the Young's modulus, or engineering elastic modulus, of the material in 
the ith principal material direction, and 1/E, is called the compliance of the material. 
Continuing along this line of thought, the off-diagonal terms of the elastic compliance, S , 
are related to the compliances, 1/E„ and the Poisson's ratios, such that Sy = -j^/E,. 
For an orthotropic material the Poisson's ratios are defined as the ratio of the strain in the 
transverse direction, j, relative to the strain in the loading directions, i, i.e. v v . = -£/£;. 

The definition of the Poisson's ratio functions for an orthotropic viscoelastic 
material are analogous to those for an orthotropic elastic material. After the manner of 
Staverman and Schwarzl, it is: 



d g i( T ) 

dr 


dr, 


(37) 


where ^(t-r) is a viscoelastic Poisson's ratio function for strains in the jth direction due to 
strains in the ith direction. Note that there is no distinction between a Poisson's ratio 
function that would be associated with a creep test or relaxation test. A proof is provided 
later demonstrating that there is no difference in the Poisson's ratio functions that result 
from a creep compliance test or a stress relaxation test. Following the initial discussions 
relative to the loading history and the nature of the convolution integral in equation (37) 
the definition may be rewritten as: 

fij = - v \> * dg i = -£■,* d *V (38) 

For the time being it will be assumed that there is no difference between Poisson's ratio 
functions resulting from stress relaxation tests and creep compliance tests. 

The relation between the off-diagonal compliance terms and the Poisson's ratio 
functions requires some development. Consider a uniaxial tensile creep compliance test in 
the first principal material direction, i.e. a x = a 0 l, and a 2 = <r 3 = 0. Substitution into 
equation (14) yields: 

£ 1 = §11 <y 0 , e 2 = §12 a of S 3 = §13 

Then, substitution of the above into equation (38) produces: 
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S 2 = S 12 a o= ~ "l2 * d£ l = “ "l2 * d t S ll J = ~ *>12 * dS„ 


e 3 — S 13 a Q — v 13 * dgj — v 13 * d[S 13 oj — o a y 13 * dS^. 

Using the definition of the tensile engineering creep compliance functions, F ; , the above 
equations yield: 

S 12 = - v n * dF! (39) 


S 13 = — p 13 * dFj. (40) 

Similarly, uniaxial tensile creep compliance tests in the second and third principal material 
directions yield, respectively: 

s 21 — ~ v 2\ * dF 2> (41) 

$23 = — ^23 * dF 2> (42) 


$31 — **31 * dF 3> 


S 32 = ~ "32 * dF 3- 


As a consequence of symmetry of the compliance tensor (matrix) and the above relations 
the i-j and j-i Poisson's ratio functions are interrelated by the following equations: 

"12 * dF i = "21 * dF 2> (45) 

"13 * dF l = "31 * dF 3> (46) 


"23 * dF 2 = "32 * dF : 


These results are analogous to those found in the development of orthotropic elasticity. 

Using matrix methods the tensile compliance tensor (matrix) may be expressed as: 
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[-Tension 


Sji s 12 S 13 

s 

= 

S 2 1 S 22 S23 

L J 


_ s 3 i s 32 S 33 _ 


= P * d F , 


1 -p n -p u 


"*T 

0 

0 “ 

-v 21 1 -P-22 

* d 

0 

F 2 

0 

1 — 

1 

J 8 

1 

K) 

1 


0 

0 

F 3 _ 


(48) 


where P is the Poisson's effect matrix and F is the tensile compliance matrix. The 
Poisson's effect matrix is a nonsymmetric matrix while the tensile compliance matrix is a 
diagonal matrix. 

Equation (20) gives the relationship between the orthotropic viscoelastic stiffness 
functions and the compliance functions and equation (22) gives the tensile or normal part 
of these relations. With equations (45) - (47), equation (22) can be rewritten as: 



Tension 


* d 


P * d F 


(49) 


By applying the identities for convolution integration the above equation may be rewritten 
as: 



* d 




(50) 


Note that the tensile compliance matrix, F, is a diagonal matrix. If we introduce the 
functions Fj(t), such that F- 1 (t) is the inverse convolution function of F ; (t), i.e. 

F;(t) * dF-^t) = 1 from convolution identity (VII), then we can write: 


^ ij * d F' 1 = 


..Tension! 


1 * d P * d F J * d F 1 , 


which, by application of convolution identity (II), becomes: 

■Tension! 


I * d F' 1 = 


J 


* d 


£pj ■ * d [f] * d F 1 


and, finally by application of convolution identity (V) and (VII), the above becomes: 


12 



(51) 


M - [ c 


.Tension 


* d 


[']• 


Note that the matrix F 1 is a diagonal matrix, while the matrices c Tension and P are full 
matrices. Further, as previously mentioned, P is a nonsymmetric matrix. The F matrix is 
eliminated from the right hand side of the above equation by the simple introduction of F 1 
only because F is a diagonal matrix. 

At this point in the elasticity development, the elasticity form of the above equation 
is simply post multiplied by the inverse of the Poisson's matrix, resulting in an expression 
for the stiffness terms as functions of the inverse compliances and the Poisson terms. The 
viscoelastic development is not as simple. While it is mathematically and analytically 
feasible to introduce a matrix II such that the convolution of II with the Poisson’s matrix, 
P, produces the identity matrix, i.e.: 

[nl • dfp] = [i], (52) 


it is computationally impractical. Allowing that the II matrix exists, by convolution of 
equation (51) with II yields: 


F 1 * dj^nj = 


pension * dj’pj * d ^ nj 


Then, after application of the convolution identities: 
^Tenskmj _ jyj . d[n], 


(53) 


in which F 1 is a diagonal matrix of inverse convolution functions of the tensile creep 
compliance functions in the principal material directions and II is a full nonsymmetric 
matrix of functions related to the Poisson's ratio functions. 

For the moment assume that the Poisson's ratio functions are known functions of 
time. To determine the elements of II requires the computation of multiple convolution 
integrals and inverse convolution functions. For proof consider the following 
development. Expand equation (52) in terms of the individual components; 


n n n I2 n 13 


1 -V n -Pi3 


'l 0 o" 

n 2 i n 22 

* d 

""21 1 -"23 

= 

0 10 

_n 3 i n 32 n 33 _ 


1 

I 

i 

to 
►— * 

1 


.0 0 1. 


13 



n„-n 12 * dv 21 -n 13 * dv 31 
n 2 i-n 22 * d^j-n^ * dv 31 
n 3 i-n 32 * dv 21 -n 33 * dv 31 


n 12 -n u * di/ 12 -n 13 * dv 32 
n 22 -n 21 * dvjj-n^ * dv 32 

n 32 -n 3 1 * dt'^-ii-j-j * dv 32 


n 13 -n„ * d»»i 3 -n 12 * d^ 

1123-112! * dP 13 -n2 2 * d^ 
n 33 -n 31 * dv 13 -n 32 * dv^ 


10 0 
0 10 
0 0 1 


Concentrating on evaluation of the n n term, expand the first row to three independent 
equations; 

n n = i + n 12 * dv 21 + n 13 * d^ 3i (i) 

n , 2 = n n * dv n + n 13 * dv 32 (ii) 

n i3 = n u * dv 13 + n i2 * dv^ (iii) 

Equations (ii) and (iii) by taking the convolution of the equations with v 13 and v n , 
respectively, thus: 

n i2 * dv 13 = (n u * dv 12 ) * dv 13 + (n i3 * dv 32 ) * dv 13 , (iv) 

n i3 * dp n = * <*" 13 ) * dp n + (F-n * d "23) * dj, i2- ( v ) 

Then the n n term is eliminated by subtracting equation (v) from equation (iv). This 
yields: 

n i2 * dv 13 — n 13 * dv 12 = 

= (IIj! * dv 12 ) * dv 13 4- (n i3 * d v 32 ) * dv X3 — (Iljj * dv, 3 ) * dv l2 ~ 

(II 12 * dv^) * dv 12 


II 12 * dp l3 — n 13 * dv 12 = (n 13 * dv 32 ) * dv 13 — (n 12 * dv^j) * dv, 2 (vi) 

Then rearranging and applying the convolution identities gives: 

n i2 * d ("l3 + *23 * ~ n i3 * d (*12 + *'32 * dv u)- ( vii ) 

Introducing the convolution inverse of (v 13 + v n * dv 12 ) it is possible to express n 12 in 
terms of II 13 and v--. This expression is: 


n i2 ~ | n 13 


* d 


v + * dv 

12 32 


l * d i * dv, L 

13 f 1 13 23 12 f 


-I 


(viii) 
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-I 

where <j v + v * dv 0 [ is treated as a single function of time such that: 


{ v,, + v ' * di» L * dJ i> + v * dp„ L 

13 23 12 f I 13 23 12 | 


-I 


= 1. 


(ix) 


By application of the convolution distribution property equation (viii) becomes: 

n 12 = n 13 * d |v 12 + * 32 * dt» 13 1 * d j v 13 + * dp 12 | . 

Equating equation (ii) and equation (x) and rearranging gives: 


(x) 


n 13 * d 


( v + v * dp L * dJ + p * dp l 
12 32 13 r 1 13 23 12 f 


-I 


- V 


32 


= n„ * d p n . 


(xi) 


Solution of (xi) for n i3 in terms of n n and p {j requires introduction of another inverse 
convolution function where 


|„i2 + , 32 * d„ 13 | * d j v 13 + ^ * d^ 12 | - , 32 j 

{"|2 + "32 * d "l3 J * 4 ",3 + "23 * d ",2 | - "32 


= 1, (xii) 


as before this function must be treated as an independent function of time. Applying the 
result from equation (xii) to equation (xiii) yields an expression for II 13 . It is: 


H13 = n n * d>{ 


12 



- 




-I 

-I 

* d 


p + p * dp 
12 32 13 

* d 

p,, + * dp,,, 

13 23 12 

— ^32 



- 

- 



• 



.(xiii) 


Then substituting the above equation into equation (x) produces: 
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I 

n i 2 = n ii * d j v n 


* d 







-I 

« 


v „ + i> * dp „ 

* d 

v „ + v * dp 

> 



12 32 13 


13 23 12 


. 

. 

. 


. 



* d 






-I 

- 



j/ + * di> 

* d 

v „ + i> * dp 

— V'xo 




12 32 13 


13 23 12 

32 



- 

- 


. 

- 

_ 



(xiv) 


Finally, substituting equations (xiii) and (xiv) into equation (i) yields an expression for n u 
in terms of p^. This expression is: 


Ml 



1 — V-, * d< 

v * d 

< 

' 

»/ + v * d v 

* d 

' 

v + * dp 


21 * 

12 

- 

> 

12 32 13 


13 23 12 


_ 

k 







-i 


* d 


— 




_ 

-I 


v + v * dv 
12 32 13 

* d 

’ ' 

v,. + » * dj> 

13 23 12 

“I 

“ "32 

- 

. 


> - 

_ 

_ 


31 



_ 







< p * d 
12 


v,. + p.. * d p„ 
12 32 13 

* d 

p t , + j> * d»> 

13 23 12 

-I 

— P 32 


► 



l i 


L 

-* 

- 



(54) 


Thus, it is possible to derive the elements of inverse convolution matrix. Computing the 
value of these elements as functions of time is practical impossibility in view of the need to 
compute three different convolution inverse functions from the known functions p--. 


2 . 5 . Stress Relaxation Testing of a Specially Orthotropic Viscoelastic Material 

The previous development focused primarily on the creep compliance formulation 
of specially orthotropic viscoelastic response. By doing so it was possible to directly relate 
the creep compliance elements, S^, to the tensile and shear creep compliance functions, Fj, 
and G u and the viscoelastic Poisson's ratio functions, u~. Additionally, it was shown that 
the shear creep compliance functions are inverse convolution functions of the shear stress 
relaxation functions and finally that, while a similar relationship exists for the tensile 
functions, the relationship is a complex combination of the Poisson's functions. Because 
the shear stress relaxation functions are directly related to the creep compliance functions 
and thus may be computed based evaluation of the convolution inverse it is not necessary 
to specifically address shear stress relaxation testing. 
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Consider a uniaxial tensile stress relaxation test for a specially orthotropic material. 
In this test a uniaxial tensile strain is applied as a step function in time in the first principal 
material directions, x t . Thus, the applied loading is: £, = £ 0 1. The boundary conditions 
for the equilibrium equations for the specimen determine the complete loading of the 
specimen. These conditions are: <r 2 = <r 3 = 0. The normal strains in the second and third 
principal material directions are not known a priori, they must be determined from the 
constitutive equation. From the previous development: 



(55) 


Then, for the stress relaxation test in the first principal material direction: 


^1 


l — 

o' 

O 

O 

0 

► z= 

C21 C-22 Qa 

0 


_ On c 32 c 33 _ 


* ds 


Co 1 

£2 >■ 
fi 3 


(56) 


To produce useful results it is necessary to rearrange the equations such that the knowns 
are grouped together on the right hand side of the expression. When this is done we find: 


L = C + Jc * dC -C * dC 


11 


23 


23 


22 


; 1 

33 [ 


-I 


* dJ 


C * d 
33 


c , * dc l2 


+ C 22 * d ( C .3 * dC ,3) - 2C ,2 * d ( C .3 * dC 23> 


(57) 


Thus, the stress measured in the stress relaxation test is related to the tensile stiffness in the 
first material direction, C n , as well as other stiffness terms. In addition, it is important to 
note that C n is not the relaxation modulus, Ej(t), in the first principal material direction. 
Referring back to equation (53) it is clear that C n is equal to F^ 1 * dn n , because F^ 1 is the 

convolution inverse of the creep compliance function it may be considered as the relaxation 
modulus function. Hence, this tensile stress relaxation test does not directly yield the stress 
relaxation modulus in the first principal material direction. Provided that other tests are 
performed such that information on C 22 , C 33 , etc. it is possible to determine the specially 
orthotropic stress relaxation functions, However, the complexity of the expression makes it 
a practical impossibility to relate the stress measured in this test to a material property 
function. 
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Viscoelastic Response 
Introduction of the isotropic behavior to the analysis of linear viscoelastic response 
generates significant simplifications in the constitutive relations. For an isotropic material 
we have the additional matrix symmetry relations: 

C = C u , for i = 1 to 3, 

G = C u , for i = 4 to 6, (i.e.: G = G 12 = G 13 = G^), 

L = Cjj, for i,j = 1 to 3, i *= j 
S = S ;i , for i = 1 to 3, (i.e.: S = F = F ; , for i = 1 to 3) 

J = Sji, for i = 4 to 6, (i.e.: J = J 12 = J 13 = J^), and 


M = Sjj, for i,j = 1 to 3, i ± j, 

where C(t) is a normal stiffness relaxation function, L(t) is a normal stiffness relaxation 
function due to Poisson’s effect, M(t) is a normal creep compliance function due to 
Poisson's effect, and G(t), J(t), and F(t) are the isotropic form of the functions previously 
introduced. With these functions and the above equalities the stiffness matrix, equation 
(17), becomes: 


Cu C 12 C 13 0 0 0 

C 12 C 22 0 0 0 

Cj 3 C* C 33 0 0 0 

0 0 0 C w 0 0 

0 0 0 0 C 55 0 

0 0 0 0 0 c 66 _ 


C L L 0 0 0 
L C L 0 0 0 
L L C 0 0 0 
0 0 0 G 0 0 
0 0 0 0 G 0 
0 0 0 0 0 G. 


(58) 


and the compliance matrix, equation (18), becomes: 


S n S I2 S 13 0 0 0 

S|2 S 22 S a 0 0 0 

S i3 S 33 0 0 0 

0 0 0 S M 0 0 

0 0 0 o s S5 0 

0 0 0 0 0 s 66 _ 


F M M 0 0 0 
M F M 0 0 0 
M M F 0 0 0 
0 0 0 J 0 0 
0 0 0 0 J 0 
_ 0 0 0 0 0 J. 


(59) 


Given the isotropic form of the compliance and stiffness matrices, the interrelation of the 
compliance and stiffness, equation (20), becomes: 
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,( 60 ) 


C L L 0 0 0 


F M M 0 0 0 


100000 

L C L 0 0 0 


M F M 0 0 0 


0 1 0 0 0 0 

L L C 0 0 0 

* A 

M M F 0 0 0 


001000 

0 0 0 G 0 0 

* a 

0 0 0 J 0 0 


000100 

0 0 0 0 G 0 


0 0 0 0 J 0 


000010 

_0 0 0 0 0 G_ 


_ 0 0 0 0 0 J. 


.0 0 0 0 0 1. 


which yields three unique equations. They are: 


C * dF + 2L * dM = 1, 

(61) 

F * dL + C * dM + L * dM = 0, 

(62) 

G * dJ = 1. 

(63) 


Equation (63) is the well known interrelation of the shear creep compliance function and 
the shear stresss relaxation function in the form given by Staverman and Schwarzl. From 
equations (61) and (62) it is possible to derive expressions for C and L in terms of F and 
M, and vice versa. These expressions are: 


C = (F + M) 


* d 


(F + 2M) 1 * 


d(F-M) , 


(64) 


L = - M * d 


(F + 2M) 1 * 


d(F - M) 


(65) 


Further, from equations (61) and (62) we find: 


(C - L) * d(F - M) = 1, 

(66) 

(C + 2L) * d(F + 2M) = 1, 

(67) 


which are easily proved from the results given in equations (64) and (65). The significance 
of equation (66) is that the inverse convolution relationship between the tensile property 
functions is held between (C - L) and (F - M) as well as between the tensile stress 
relaxation modulus function, E, and the tensile creep compliance function, F. 

Understanding of the above isotropic equations is expanded by the introduction of 
the viscoelastic Poisson's ratio function. The orthotropic form of the Poisson's ratio 
function is introduced in equation (37), the isotropic form is simply: 



dfii(T) 

KM— |^-dr. 


-00 


(67) 
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or in convolution notation: 


£j = - v * d£ ; = - gj * dv. (68) 

As a result of this definition and the results from the orthotropic development the off 
diagonal terms of the compliance matrix are written as: 


Sy = M = - v * dF, for i,j = 1 to 3, i j. 


(69) 


With this identity and following some manipulation, in which convolution identity (Vm) is 
of particular importance, equations (64) and (65) become: 


C = F 1 * d 


-i 


(1 — p) * d (1 + v) * d(l — 2v) 


-i 


(70) 


and 



v * d 


(1 + v) 1 * d(l- 2v) 


(71) 


Note that if the material is elastic, i.e. i>(t) = v t , a constant the above equations reduce 

E(1-jO 

to the stiffness values found in the generalized Hooke's law; C = — — ~ — : , — - and L = 

(l-2i/ e )(l+p e ) 


(l-2v e )(l+p e y 


where E is the Young's or elastic modulus. 


Thus, providing reassurance 


of the viscoelastic results. 

Returning again to the discussion of a uniaxial tensile stress relaxation test, consider 
the results produced from such a test for an isotropic material. Using the stiffnesses, C 
and L, and an applied strain history: e x = £ 0 1, cr 2 = cr 3 = 0, (e 2 , £ 3 £ 0), we find: 


= 


C - 2(L * dL) * d(C + L) 


(72) 


and 


£ 2 = £ 3 = - L * d(C + L) 1 £ 0 (73) 

or in terms of the Poisson's ratio function: 

= F" 1 £ 0 ( 74 ) 


and 



&2 £3 


(75) 


= —v e 


o' 


The consequence of which is the function found by dividing the measured stress, o x , by the 
applied strain magnitude, e 0 , is the tensile stress relaxation modulus of the material. 
Further, with the introduction of the tensile stress relaxation function, E(t), it is clear that 

E(t) is equal to F (t), and therefore: 

F * dE = 1, (76) 


which is a classical result and is stated without proof by Staverman and Schwarzl. Using 
the tensile stress relaxation function, E(t), the expressions for the tensile stiffness 
components, C and L, equations (70) and (71), can be rewritten as: 


C = E * dj 


(1 — v) * d 


(1 + v) 1 * d(l— 2v) 1 


(77) 


and 


L = E * d 


-1 


* d (1 + v) * d(l- 2v) 


-1 


(78) 


These expressions allow direct computation of the two tensile stiffness functions from the 
stress relaxation function and Poisson's ratio function. It is not apparent from equation 
(74) or (75) would reduce to a simpler form without the introduction of the Poisson's 
ration function. 


2.7. Poisson's Ratio Function in Isotropic Viscoelastic Materials 

The definition of the Poisson's ratio function given for specially orthotropic 
material in equation (37) and for an isotropic material in equation (67) do not place any 
restrictions on the applied normal strain. Staverman and Schwarzl state that there are a 
number of possible definitions of the Poisson's ratio function for a viscoelastic material. 
They further state that: "two possibilities are: 1. The ratio of the lateral contraction to 
longitudinal elongation in a creep experiment. 2. The ratio of lateral contraction to 
longitudinal elongation in a stress relaxation experiment. For a purely elastic material both 
definitions are identical. For a viscoelastic material the second possibility has bee chosen 
rather arbitrarily." A proof that in viscoelasticity, as in elasticity, there is only one 
Poisson's ratio function is presented in this section. 

For this proof, the discussion is restricted to isotropic viscoelastic materials. 
Recalling the definition of the Poisson's ratio function for an isotropic material from 
equation (67) write: 


21 



£j(t) = - 


dfii(r) 

vS- T ) ~ — dr 


dr 


-oo 


and 


£j(t) = - 


d£;(r) 

"rC-*) ~n~ dr, 


dr 


-00 


(79) 


(81) 


where ^ c (t) is the Poisson's ratio function that results from a creep compliance test and 
i> R (t) is the Poisson's ratio function that results from a stress relaxation test. Again, for 
convenience use the convolution notation for the above equations, thus: 

£j = - * d£ ; = - £; * dv c . 


and 


£j = — * dfij = — £ ; * dp R . 

The discussions presented previously for creep compliance testing of a specially orthotropic 
viscoelastic material and isotropic viscoelastic material are directly applicable to this 
discussion and the previous conclusions may be used directly with the generalization to 
isotropic materials. Therefore, using equation (70): 

Sjj = M = — v c * dF, for i,j = 1 to 3, i j. (82) 

Define that the off diagonal components of the tensile stiffness matrix, for isotropic 
materials this is L, are determined from a stress relaxation test; thus, L = L(v R ). In the 
previous development for a uniaxial tension stress relaxation test it was found that: 

-i o x 

F = E — (83) 

fc 0 


and 


j> r = L * d(C + L) \ (84) 

where j> r is used here to indicate the result is developed from a stress relaxation test. 
Rearranging equation (84) and solving for L in terms of and C, and subsequently 
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substituting this expression into equation (73) and using (83) results in an expression for C 
in terms of v R and E, derived solely from a tensile stress relaxation test. The result is: 


C = E * d 


[< 


a - -r) * dla + XRf' » d(l- 2 (' r )~ i 


(85) 


When the above result is substituted back into the intermediate expression for L an 
expression for L in terms of v R and E, derived solely from a tensile stress relaxation test, is 
obtained. This expression is: 


L = E * d 




"r * d I (1 + *>r) 1 * d(l- ^r) 1 


( 86 ) 


Equation (66) is a general expression that interrelates the tensile viscoelastic material 
property functions. It is: 


(C - L) * d(F — M) = 1. 

Using equation (82) in rewriting the above relation yields: 


( 66 ) 


C - L 


* d 


F * d 


1 + V r 


= 1 . 


The results from the tensile stress relaxation test, equations (85) and (86) are now 
substituted into the above equation. This yields: 


= 1 . 


From the interrelation of tensile stress relaxation modulus and tensile creep compliance, 
equation (77), E and F are eliminated from the above equation. Therefore: 


E * d 

fl+,J _I * d (l - 2j/ r ] 1 

V J l J 


* d| F * d 
l 

1 + v c 




4 





1 + Vr 


= + "R > 


which simplifies to: 

v c = J'r. (87) 

Thus, the Poisson's ratio function found from a creep compliance test in tension, v c , is the 
same as the Poisson's ratio function found from a tensile stress relaxation test, y R . The 
identical result, equation (87), is obtained from the general expression in equation (67). 
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To conclude then, there is only one Poisson’s ratio function in linear isotropic 
materials. It may be measured by either a tensile stress relaxation test, or by a tensile 
creep compliance test. 
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3. An Alternate Form of the Deformation Invariants 

The general purpose finite element code Abaqus uses an alternate form of the 
deformation invariants in the constitutive equations implemented in the code. Dr. S. Peng 
has proposed using a similar alternate formulation for the constitutive equation that he is 
developing. This section addresses these alternate forms and there relationship to each 


other and the hyper-elastic constitutive equation. 

Define the quantity J as: 

J = det( F), (1) 

where F is the deformation gradient. Introduce a new measure of the deformation, 4>. 
Where the new measure is defined as: 

* = J~ 1/3 F. (2) 

The Left Cauchy-Green Deformation Tensor, B, is defined as: 

B = FF‘. (3) 

From this definition a new deformation tensor, B, is defined as: 

B = ** T . (4) 

Therefore: 

B = J m B, (5) 

or 

B = J m B. (6) 

The Cayley-Hamilton Theorem for the tensor B states: 

B 3 - I,B 2 + 1^ - y = 0, (7) 

where I, are defined as the invariants of the tensor B. The compatible definition of the 
invariants is: 

I, = 'KB), (8) 

t, = J( I, ’ - (KB’)) , (9) 

t, = dam. (10) 
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Note that the new deformation measure, B, is a symmetric tensor and therefore must also 
satisfy the Cayley-Hamilton Theorem. Therefore: 

B J -r 1 B 2 + r 2 B-r 3 i = o, (ii) 


where T. are defined as the invariants of the tensor B. Similarly, the new invariants are: 

r, = tri B), (12) 

r, = \ /rf - 

r, = det( B). 

Substitution of the result from equation (6) into equation (12) with equation (8) yields: 

T, = rKJ'^B) = J M /r(B) = T 2 *!,. (15) 


(13) 

(14) 


Similarly, for r 2 : 


r, = j Jr 3 !; - r“V(B) ! J, 


r = — 

2 


fi - (KB 2 )} = 1"% 

and for T 3 : 

r 3 = detii^B) = (Tydeti B) = (J 2/3 ) 3 I 3 . 

Recall equations (1), (3), and (10), they are combined to give: 
L, = det(B) = det(FF J ) = det(F)det(F J ) = det(F)det(F) 
I 3 = J 2 . 

As a result, then: 

r 3 = (O 3 J 2 = 1. 


(16) 


(17) 


(18) 


(19) 


As a proof that the transformation of the deformation measure is consistent consider 
equation (7), the Cayley-Hamilton Theorem for the original deformation tensor: 
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B 3 - I,B 3 + I^ 3 - y = 0, 


substitution of the definition of the new deformation measure yields: 
(j w B) 3 - i.cTb) 2 + y w B - y = 0, 

B 3 - Iir w B 2 + I 2 r"B - I 3 J^I = 0 . 


Then, from equations (15), (16), and (17): 


b 3 - r,B 2 + r 2 B - r 3 i = o, 


which is equation (11) the Cayley-Hamilton Theorem for the new deformation measure. 
Therefore, the definitions of the new deformation measure and the invariants are 
consistent. 

The constitutive equation for a compressible material is given as: 


a = 2r 3 ,/2 {l 3 W 3 I + (W, + I.WJB - W 2 B 2 }, 


( 20 ) 


where a is the stress, and W, is the partial derivative of the strain energy density function, 
W, with respect to the ith invariant of the deformation, I,, i.e.: 


awq.,1,,13 

31, ’ 

awq,,yy 

a, ’ 


(21) 

( 22 ) 




(23) 


Now introduce the new deformation measure into the constitutive equation: 
a = 2 t n jW + (W, + I.W^J^B - W^J^B) 2 j , 

a = 2 ll a J 4 * |l 3 W 3 J' 4 ' 3 I + (W, + yV^B - W 2 B 2 j, 

and with the new invariants: 

a = 2(j 6/3 r 3 )- , ' 2 r|(j 6/3 r 3 )w r r 3 i + [w, + (j M r,)wjj M B - w 2 b 2 J, 

a = 2I T 3 ' 2 J 1/3 f(f 3 r 3 )W 3 I + [T^W, + (rr,)W 2 ] B - W 2 B 2 l. (24) 
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Note that W. are the partial derivatives of the strain energy density function with respect to 
the original invariants. To complete the transformation process introduce Q, the strain 
energy density function expressed as a function of the new variables. Thus: 


wo,,!,,:,) = w(j“r„i“r ! ,rr 1 ) = fi(r„r„r„j). 

Therefore, the W ; become: 

_ awq, _ 3fl(r„r„r„j) ar, d n dJ 
Wi ai, ar, ai, + aj ai,’ 


(25) 


w, = 


an(r„r 2 ,r 3 ,j) 
ar, 


f = J Q„ 


(26) 


_ dwg^D _ a0(r„r 2 ,r 3 ,j) ar 2 afl(r„r 2 ,r 3 ,j) aj 
Wz ar 2 ai, + aj ay 


w. 


= 3tl(r ':5’ r, ' J) r = r n. 


ar 


(27) 


aw([„i,,ij aO(r„r„r„j) ar, ao(r„r„r„j) ar, ai)(r„r„r„j) ar, 
W) a, ar, ai, ar, ai, ar, ai, 

, ati(r„r„r„j) ai 


aj ai,' 


ao(r„r„r„j) r j-' . 
w >~ ai 1~ 2 


(28) 


Substitution of equations (26), (27), and (28) into equation (24) yields: 


a = 2YW 3 


a — 2F 3 i/2 J 1/3 J 


(j^rjy n,i + [r“j“o, + (j 4/3 r,)j 4/3 oj b - rt 2 B 2 j, 

r l/3 r n 

yr^i + jo, + j^r.nJ b - j 4/3 q 2 e 


and since T 3 = 1, 

a = j‘* Jr m n ; i + 2 j 8/3 [ + r,n 2 ] b - 2 j 4/3 q 2 b 2 J, 

a = 0,1 + 2[T 1/3 n, + J 3 r,n 2 ] B - 2f%B\ (29) 
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» - * 


* 


with the strain energy density function formulated as: 0 = Q(T,,r 2 ,J). It is not 
immediately apparent that the expression for the stress in equation (29) is significantly 
simpler than that given in equation (20). Both forms of the strain energy density functions, 
W and 0, are dependent on three variables. Thus, the problem of identifying the strain 
energy density function is not greatly simplified. Equation (29) does begin to look like the 
expression for the stress in an incompressible material. In fact for an incompressible 
material, L, = 1 = J. Then equation (29) becomes: 

<r = 0,1 + 2 [0, + r,Q 2 J B - 20 2 B\ (30) 

and 

a = W 3 I + 2 jw, + r,W 2 j B - 2WjB 1 2 , 

a = W,I + 2 jw, + I,W 2 ] B - 2W 2 B\ 

a = W 3 I + 2 |w, + I,W 2 J B - 2W 2 B 2 . 

For an incompressible material the constitutive equation is: 
a = -pi + 2 [w, + I,W 2 1 B - 2W 2 B\ 


Therefore, for I, = J = 1, Wj^,,!,,!) = 0 3 (r,,r 2 ,l) = -p. Where p is an arbitrary 
hydrostatic pressure that can only be determined by application of the equations of motion 
and the boundary conditions. 


1 Staverman, A. J. and Schwarzl, F., "Linear Deformation Behavior of High Polymers," in Die 

Physik Der Hochpolymeren, Vol. 4, ed. by H. A. Stuart, Springer-Verlag, Berlin, 1956. 
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